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FOREWORD 


This  final  report  is  a  documentation  of  the  work  done  under 
contract  F44620-67-C-0080,  Investigation  of  PKP  Seismic  Waves.  Details 
of  the  project  history,  data  gathering  operations,  and  results  of 
preliminary  examination  of  the  data  were  presented  in  an  Annual  Report 
issued  March  15,  1969,  and  are  only  summarized  in  this  report.  Project 
results  and  conclusions  are  presented  here. 

Stanford  Research  Institute’s  participants  in  the  project  were  as 
follows:  C.  W.  Smith,  Physicist;  W.  F.  Isherwood,  Geophysicist; 

D.  E.  Wolf,  Mathematician;  and  J.  D.  Eisler,  Senior  Geophysicist. 

Initial  project  supervisor  was  W.  H.  Westphal;  interim  supervisor  was 
R.  A.  Schmidt,  and  supervisor  for  the  conclusion  of  the  project  is 
D.R.  Grine,  Manager,  Geophysics  Program.  Major  D.  D  Young,  AFOSR, 
acted  ao  initial  Project  Monitor  and  also  assisted  with  the  installation 
of  the  seismic  array  at  Byrd  Station. 

Bruce  A.  Bolt,  Professor  of  Geophysics,  University  of  California, 
Berkeley,  acted  in  a  consulting  capacity. 
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INTRODUCTION 


Stanford  Research  Institute  (SRI),  with  support  from  the  Advanced  Research 

Projects  Agency  (ARPA  and  the  National  Science  Foundation  (NSF) ,  installed 

a  four-element  seismic  array  (BYA)  at  Byrd  Station,  Antarctica,  beginning 

operation  January  18,  1968.  Short  period  vertical  component  seismometers 

were  selected  in  anticipation  of  high  magnification  operation.  Figures  1 

and  2  show  the  plan  of  the  BYA  array  and  a  block  diagram  of  the  system. 

The  March  15,  1969,  Annual  Report  (Isherwood,  1969)  describes  in  detail 

the  system  and  the  data  gathering  operations.  The  site  was  selected  to 

study  core  phase  arrivals  because  PKP  waves  have  a  principal  focus  at  an 

angular  distance,  A,  of  approximately  142°  from  the  source,  and  Byrd 

o 

Station  is  close  to  142  from  several  of  the  great  seismic  zones  in  the 
Northern  Hemisphere,  as  shown  in  Figure  3. 
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FIGURE  1  PLAN  OF  BYA  SEISMIC  ARRAY 
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FIGURE  2  BLOCK  DIAGRAM — TYPICAL  CHANNEL 


FIGURE  3  142°  GEOCENTRIC  OISTANCE  FROM  BYRD 


Helicorder  and  magnetic  tape  ecords  for  the  period  January  18,  1968, 
through  November  24,  1969,  including  approximately  1000  core  phase  arrivals, 
were  analyzed.  The  objectives  of  the  research  contracted  between  ARPA  and 
SRI  were  as  follows: 

1.  To  establish  and  prepare  specifications  for  a  high-quality 
seismograph  station  capable  of  recording  PKP  seismic  waves 
to  operate  in  Antarctica  during  the  year  of  1968, 

2.  To  study  the  characteristics  of  the  ambient  seismic  noise 
and  relate  them  to  oceanic  and  continental  storms. 

3.  To  apply  appropriate  data  processing  techniques  aimed  at 
enhancement  and  identification  of  PKP  waves  in  the  presence 
of  ambient  seismic  signals  and  noise. 

4.  To  analyze  and  interpret  PKP  wave  data  in  terms  of  their 
source  characteristics  and  the  internal  constitution  of 
the  earth. 
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Since  the  above  Interpretations  depend  on  accurate  travel  time  prediction, 
BYA-observed  travel  times  were  compared  with  published  tables  as  an 
inseparable  objective. 

Initial  evaluation  was  based  on  the  origin  times,  hypocontors,  and 
magnitudes  published  by  the  U.  S.  Coast  and  Geodetic  Survey  (USCfcGS). 
Events  listed  herein  are  referred  to  by  the  USCItGS  origin  times. 

Nomenclature  will  follow  that  of  Jeffreys  and  Bullen  (1967)  and 
Bolt  (1968).  Figure  4  presents  the  PKP  travel  time  curve  for  reference 
throughout  the  text. 


RESULTS 

Seismic  Noise 

Although  background  noise  level  varied  with  time,  Byrd  Station  was 
found  to  be  a  fairly  quiet  seismic  site,  especially  within  the  band  of 
PKP  frequencies.  Hellcorder  records  were  generally  made  with  a  system 
magnification  near  100  K.  Typically,  higher  noise  levels  were  observed 
for  periods  of  one  to  two  weeks  Interspersed  with  longer  periods  of  low 
noise  levels. 

The  seismometers  were  buried  at  a  depth  of  ~  60  m,  providing  adequate 
insulation  from  direct  vibration  by  local  winds.  Inspection  of  the  records 
showed  no  correlation  between  the  seismic  noise  level  and  local  storms; 
many  severe  local  storm  occurred  during  periods  of  low  noise.  The  months 
of  March  and  April  1968  present  an  example:  During  March,  the  records 
were  generally  noisy  with  the  most  noise  from  about  the  10th  to  19th. 

Wind  records  for  this  period  show  a  two-day  storm  on  the  13th  and  14th 
but  show  lower  average  winds  and  peak  gusts  than  later  in  the  month. 

April,  which  set  a  record  monthly  average  wind  speed  of  27.1  mph,  had 
relatively  low  noise  levels. 

Correlation  of  noise  levels  with  oceanic  storms  was  prevented  by 
lack  of  meteorological  data  for  the  oceans  near  the  Byrd  Land  coast. 
According  to  the  Antarctic  Support  Activities  (OIC,  Detachment  Charlie, 
personal  communication),  there  were  a  maximum  of  five  ships  in  the  entire 
area  off  the  coast  in  the  calendar  year  1968,  providing  almost  no 
meteorological  data.  State-of-sea  determination  from  satellite  pictures 
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FIGURE  4  BYA  TRAVEL  TIMES 


is,  as  yet,  an  unperfected  art.  Consequently,  the  surface  charts  available 
have  a  noticeable  lack  of  data  for  the  required  ocean  area  (Figure  5). 

Tests  were  run  to  determine  the  frequency  content  of  the  ambient 
noise  and  the  PKP  signals.  Initial  frequency  measurements  were  made  by 
scaling  the  helicorder  records.  Later,  analog  filtering  of  the  magnetic 
tapes  were  employed  and  finally  spectral  analysis  was  carried  out  by 
digital  computer.  Figure  6  is  a  plot  showing  the  frequency  spectrum  of 
an  average  noise  sample  compared  with  a  typical  sample  of  PKP  signal. 

This  clearly  shows  that  the  majority  of  signal  energy  is  contained  between 
0.4  Hz  (2.5  sec  period)  and  2  Hz  (0.5  sec  period).  The  background  noise 
lies  mostly  at  much  higher  frequencies  and  slightly  lower  frequencies  than 
this  signal  range.  Figure  7  shows  an  extreme  case  before  and  after  digital 
filtering  with  pass  band  of  0.4  to  2.0  Hz.  The  large  amplitude  noise  near 
5  Hz  on  Figure  7  is  encountered  when  snow  cats  are  in  operation  in  the 
station  area.  Being  well  out  of  the  signal  band,  it  is  the  easy  target 
of  additional  filtering.  Natural  ground  motion  presents  more  of  a  problem. 
The  greatest  natural  interference  arises  from  longer  period  noise,  which 
rises  to  ~  30  mp.  at  0.3  Hz  (as  measured  from  the  helicorder)  and  may  be 
as  high  as  ~  10  ng  at  0.5  Hz.  Signals  from  1  to  2  Hz  have  least  competit¬ 
ion  with  background  noise,  which  is  generally  <  1.5  np,  of  ground  motion. 

Detection  Threshold 

PKP  arrivals  were  observed  at  Byrd  Station  for  tne  range  of  A  =  110 
o 

to  173  .  The  amount  of  energy  transmitted,  however,  varies  considerably 
over  this  range  as  pointed  out  by  Shahidi  (1968).  At  distances  where 
more  than  one  arrival  are  seen,  the  relative  energy  from  different  wave 
paths  also  varies  with  A.  This  is  seen  both  in  the  amplitudes  of  arrivals 
and  the  quality  of  the  first  motion  seen.  The  variation  is  shown  in 
Figure  8,  a  composite  of  PKP  arrivals.  Table  1  lists  the  events  used 
in  Figure  8.  Although  an  attempt  was  made  to  use  events  of  near  the 
same  magnitude  and  use  records  of  near  the  same  gain,  the  relative 
amplitudes  between  different  traces  means  little  for  such  a  small  sample. 
Events  with  A’s  <  140°  show  early  arrivals,  corresponding  to  the  GH 
branch  of  the  travel  time  curve.  These  ray-paths  transmit  so  little 
energy  that  arrivals  are  seen  only  from  large  events  and  even  then  onsets 
observed  were  emergent. 
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FIGURE  6  SPECTRAL  CONTENT  OF  NOISE  AND  PKP  SIGNAL 
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FIGURE  8  PKP  COMPOSITE 


FIGURE  8  PKP  COMPOSITE  concluded 


Table  I 

EVENTS  SHOWN  ON  FIGURE  8,  PKP  COMPOSITE 


A° 

Region 

Date 

mb 

Depth  (km) 

173.2 

Novaya  Zemlya 

Nov.  7, 

1968 

6.1 

0 

166.1 

Laptev  Sea 

Apr.  7, 

1969 

5.5 

33 

153.4 

Iceland 

Dec.  5, 

1968 

5.5 

5 

149.1 

E.  Kazakh 

Sep.  5, 

1968 

5.5 

0 

145.1 

S.W.  Russia 

Sep . 26 , 

1969 

5.6 

0 

142.3 

E.  Kamchatka 

Jan. 22, 

1969 

5.5 

33 

140.2 

Turkey 

Sep.  3, 

1968 

5.7 

5 

138.9 

Kamchatka 

Jun. 15, 

1968 

5.4 

39  D 

137.5 

Kurile  Islands 

Oct . 24, 

1968 

5. 5 

35 

135.8 

Hindu  Kush 

Dec  .19, 

1968 

5.4 

151 

133.0 

Hokkaido,  Japan 

Aug.  7, 

1968 

5.6 

54  D 

129.7 

E.  of  Honshu,  Jap 

an  Jun. 12, 

1968 

5.7 

36  D 

125.7 

Shikoku,  Japan 

Dec . 11 , 

1968 

5.4 

32  D 

120.7 

Ryukyu  Islands 

Nov. 12, 

1968 

5.8 

48  D 

117.1 

Taiwan 

Feb. 26, 

1968 

6.0 

24 

114.4 

Philippines 

Nov. 10, 

1968 

5.2 

33 

110.8 

Luzon,  Philippine 
_ 1 

s  Aug .  3 , 

1968 

5.9 

37 

The  PKP  detection  threshold  at  Byrd  Station  was  determined  as  a 
function  of  A  by  searching  the  helicorder  records  for  about  500  USC&GS 
listed  events.  The  curves  on  Figure  9  show  the  detection  threshold  and 
its  variation  with  A  from  these  unfiltered  helicorder  records.  (  Not  all 
of  the  arrivals  identified  were  recognized  by  the  on-s^te  operator.)  As 
would  be  expected  there  is  some  scatter  about  the  lines  given.  Beyond 
150  ,  the  number  of  events  listed  by  the  USC&GS  were  fewer;  consequently, 
this  portion  of  the  figure  is  poorly  determined.  Beyond  140°,  the  arrivals 
from  different  branches  frequently  interfere  with  each  other.  A  separate 
curve  for  EF  arrivals  is  shown  where  they  are  distinct.  This  threshold 
could  no  doubt  be  lowered  by  providing  sharper  cutoff  analog  filters. 
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The  threshold  determined  from  the  helicorder  records  is  not  the 
optimum  available  from  the  recorded  data.  Work  is  needed  to  optimize 
the  data  and  assess  the  realizeable  detection  threshold.  The  magnetic 
tape  data  may  then  be  machine  processed  to  lower  the  threshold. 

Data  Processing  Techniques 

Although  the  present  analog  system  is  highly  filtered  (Figure  10 
shows  the  typical  frequency  response),  the  value  of  further  data  pro¬ 
cessing  is  clearly  shown  in  Figure  7.  For  computer  processing,  the 
records  from  142  events  were  selected  and  digitized  by  the  Vela 
Digitizing  Center  at  50  samples  per  sec.  The  first  enhancement  of  PKP 
arrivals  was  by  digital  filtering  with  a  band  pass  of  0.4  to  2.0  Hz,  the 
results  of  which  were  shown  previously.  These  cutoff  frequencies  were 
chosen  by  inspection  of  the  seismograms.  No  doubt  even  greater  enhance¬ 
ment  can  be  accomplished  with  optimum  filtering — which  has  not  been 
attempted  at  this  stage. 

Summations  of  the  signals  from  the  four  seismometers  were  formed, 
as  well  as  "lagged  sums"  as  a  test  of  velocity  filtering  feasibility. 

Results  indicated  that  a  straight  sum  of  the  signals  was  helpful  in 
correlation  of  significant  arrivals,  but  the  BYA  array  was  too  small 
to  show  additional  improvement  by  velocity  filtering.  Uncertainty  of 
local  geologic  influence  and  relative  timing  between  channels  added  to 
the  difficulty  in  interpreting  values  of  apparent  velocity  observed 
across  the  array.  (The  average  of  arrival  times  from  the  playback  of 
100  events  shows  that  arrivals  on  the  center  seismometer  are  seen  0.1 
sec  earlier  than  on  the  surrounding  three  seismometers.  Geologic  and 
elevation  uncertainties  and  measured  relative  phase  delays  in  the  electronics 
cannot  account  for  this  difference.  At  these  slow  tape  speeds,  however, 
relative  timing  is  quite  sensitive  to  precise  tape  head  spacing  and  tape 
stability.  The  center  seismometer  was  recorded  on  a  different  head  from 
that  used  on  the  other  seismometers.  This  should  be  avoided  if  a  tenth 
of  a  second  reliability  in  timing  is  expected.) 
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Travel  Times 

Existing  travel-time  tables  formed  a  starting  point  for  analysis 
of  BYA  observations.  Jeffreys  and  Bullen  (J-B) ,  1940,  published  tables 
of  predicted  PKP  travel  times  based  on  a  two-layered  core  of  the  earth. 

This  included  predicted  times  for  paths  corresponding  to  the  AB,  BC,  and 
DEF  branches  of  the  travel  time  curve.  As  more  observational  data  became 
available,  additional  arrivals  (precursors)  were  noted  and  slight  deviations 
from  the  J-B  tables  were  seen.  In  1968,  new  tables  by  Bolt  were  included 
in  a  new  set  of  seismological  tables  for  P  phases  (Herrin  et  al.,  1968). 
These  included  times  for  a  new  branch,  GH;  a  consequence  of  the  inter¬ 
mediate  core  in  a  tripartite  core  model.  New  travel  times  for  the  BC 
branch  were  not  included.  These  1968  tables  (corrected  2.0  sec  for  use 
with  J-B  P  tables)  were  used  in  this  study. 

Appendix  A  lists  good  quality  PKP  arrivals  at  Byrd  used  in  this 
analysis.  Observed  travel  times  were  corrected  for  the  earth's  ellipticity 
according  to  tables  by  Jeffreys  and  Bullen  (1967).  Depth  corrections 
were  taken  from  Bolt  (1968). 

Herrin  and  Taggart  (1968)  give  a  station  correction  for  station  BYR 
(USC&GS  station  located  within  100  m  of  the  center  BYA  seismometer)  of 
0.24  sec  +  0.24  sln(Z  +266°),  based  on  90  observations.  Simple  calculation 
of  the  expected  delaying  effect  of  the  seismometer's  location  on  2160  m 
of  ice  and  at  an  altitude  of  1500  m  above  sea  level  yields  a  station 
correction  of  0.43  sec.  An  observational  check  of  82  teleseismic 
P  (75°  <  A  <  95°)  arrivals  at  BYA  during  1968  showed  an  average  residual 
of  0.8  sec  (compared  to  J-B  tables).  No  significant  azimuthal  dependence 
was  seen. 

Due  to  these  discrepancies,  no  station  correction  is  introduced. 

Future  evaluation  of  the  data  may  justify  further  consideration  of  this 
factor. 

o 

Between  about  142  to  152  ,  relatively  high  amplitude  arrivals  are 
identified.  These  form  a  consistant  set  with  sharp  onsets  (see  Fig.  8) 
but  do  not  fall  on  any  of  the  branches  listed  in  the  1968  tables.  On 
the  basis  of  the  recent  literature,  these  will  be  referred  to  as  BC 
arrivals  (e.g.,  Sacks  et  al.,  1970,  and  Engdahl  and  Felix,  1970). 
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These  are  plotted  separately  on  Figure  11.  Each  point  plotted  represents 
the  corrected  travel  time  from  a  single  event.  Circles  are  earthquakes; 
triangles  are  suspected  blasts  (0  km  depth  by  USC&GS) .  From  a  total  of 
115  arrivals,  a  curve  was  computed  by  at  least  squares  fit.  The  first- 
order  fit: 

t  =  19m  31.07726  sec  +  3.00003  (A  -  142°)  sec 

standard  error  =  0.746  sec 

is  shown  in  the  figure.  Second-  and  third-order  curves  were  also 
calculated  and  ure  given  below: 

t  =  19m  30.749  sec  +  3.2447  (A  -  142°)  sec  -  0.0282003  (A  -  142°)2sec 

standard  error  =  0.726  sec 

t  =  19m  30.885  sec  +  2.988  (A  -  142°)  sec  +  0.04730  (A  -  142°)2sec 

-  0.00573871  (A  -  142°)3  sec 

standard  error  =  0.724  sec 

Because  of  the  small  difference  in  standard  error,  the  second-  and  third- 

order  fits  are  not  justified  at  this  time. 

Although  some  investigators  have  extended  the  BC  branch  to  155°  or 

farther,  the  sparse  BYA  data  do  not  confirm  extension  beyond  152°.  The 

USC&GS  shows  few  earthquakes  at  these  distances;  also,  the  detection 

threshold  appears  to  be  higher  for  this  range.  As  seen  in  Figure  8, 

the  BC  arrival  has,  at  the  least,  lost  its  preeminence. 

Arrivals  identified  as  falling  on  the  AB  branch  were  observed  at 

A  >  145°,  From  145°  <  A  <  152°,  however,  the  AB  arrivals  follow  closely 

o  o 

the  larger  BC  signals.  From  153  <  A  <  180  ,  AB  arrivals  are  the  largest 

PKP  amplitudes  recorded  at  BYA.  Since  the  detection  threshold  of  AB  was 
lower  than  that  of  EF  in  this  range,  enough  AB  arrivals  were  recorded  to 
confirm  the  trend  and  general  shape  of  the  published  curve.  BYA  obser¬ 
vations,  however,  were  usually  later  than  predicted  by  several  seconds. 

For  instance,  the  well-recorded  events  originating  at  Novaya  Zemlya  give 
AB  travel  times  of  1.3  and  1.5  sec  greater  than  predicted.  Combined  with 
the  EF  arrival  times  previously  mentioned,  this  gives  DF-AB  intervals  of 
3.6  to  3.9  sec  greater  than  predicted  and  show  that  neither  can  be  adjusted 
for  a  different  depth  of  hypocenter  without  increasing  the  deviation  of 
the  other. 
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FIGURE  11  BYA  TRAVEL  TIMES  FOR  BC  BRANCH 
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Of  additional  interest  on  these  events  at  A  =  173°  is  seismic 
energy  starting  at  least  25  sec  earlier  than  the  AB  arrival  (Figure  12). 
This  does  not  appear  to  be  described  in  the  literature;  however,  Bolt 
(personal  communication)  reports  indications  of  this  energy  from  other 
events  at  similar  A's,  and  Ergin  (1967)  shows  numerous  arrivals  between 
the  DF  and  AB  branches. 

Source  Characteristics 

Part  of  the  study  of  source  characteristics  has  shown  that  core 

phases  offer  some  promise  in  discriminating  between  earthquakes  and 

explosions.  PKP  and  PKIKP  can  be  used  in  most  of  the  same  ways  as 

teleseismic  P.  As  pointed  out  by  Bolt  (1968),  they  have  an  advantage 

in  that  PKIKP  has  a  sharp  onset  beyond  A  =  160°.  PKP  has  a  large 

o  o 

amplitude  with  a  sharp  onset  for  142  <  A  <  152  on  the  BC  branch.  Bolt 

also  notes  that  the  contrast  in  dt/dA  between  the  AB  (denoted  by  PKP  ) 

and  EF  branches  for  A  >  160  produces  a  PKP  -  PKIKP  internal  more 

sensitive  to  changes  in  A  than  the  S-P  interval  near  A  =  50°.  An  error 

of  0.1°  in  A  should  be  detectable. 

Examples  are  shown  for  earthquakes  and  Russian  nuclear  shots  of 

first  motion,  determination  of  depth  from  BC  branch  travel  times,  and 

the  sharpness  of  the  PKP2  -  PKIKP  interval. 

The  maximum  sensitivity--or  lowest  detection  threshold--appears  to 

fall  on  the  BC  branch  at  142°  <  A  <  152°.  Within  this  range, numerous 

events  of  USC&GS  magnitude  m^  =  4.5  and  occasionally  less  are  seen  even 

o 

on  unenhanced  records.  Below  142  ,  there  appears  a  gradual  dropping  off 

o 

in  amplitude  and  quality  to  ~  110  .  For  example: 

117°  Nevada  Test  Site  (Boxcar)  m.  =  6.3  Yielded  15mp.  Poor  onset 

D 

o 

147.6  North-Atl antic  earthquake  m^  =  4.7  Yielded  40iHi  Good  onset 

As  for  P,  the  first  motions  of  PKP  events  are  possible  distinguishing 
characteristics  between  explosions  and  natural  earthquakes.  Figure  13 
shows  onsets  of  arrivals  of  the  BC  branch.  Trace  1  shows  the  common  case 
in  which  explosions  show  a  clearly  compressional  first  motion.  However, 
Trace  2  shows  that  these  are  not  always  clear-cut,  especially  if  only  one 
seismometer  trace  is  examined.  (In  this  case,  comparison  with  traces  from 
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FIGURE  12  DF  -  AB  INTERVAL 


COMPRESSION 


#1 

E.  KAZAKH  BLAST 

mb  “  6.4 

A  *  148.9 

#2 

E.  KAZAKH  BLAST 

mb  -  5.0 

A  -  149.1* 

#3 

N.  ATLANTIC  QUAKE 

rrib  -  4.7 

A  -  147.6* 

TA-6495-19R 


FIGURE  13  FIRST  MOTIONS 
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the  other  recording  seismometers  shows  the  true  compressional  onset.) 
Although  none  of  the  BC  arrivals  from  blasts  examined  showed  clear  dilation, 
such  as  shown  on  Trace  3  from  a  North  Atlantic  earthquake,  onsets  from 
events  of  m^  <  5.0  were  unclear  from  preliminary  examinations. 

The  sharp  onsets  observed  on  the  BC  branch  arrivals  lead  to  another 
distinguishing  criterion — depth  of  hypocenter.  Since  PKP  waves  travel 
along  paths  nearly  normal  to  the  earth's  surface,  differences  in  arrival 
times  vary  more  with  hypocenter  depth  than  for  other  wave  paths.  For 
example,  if  the  true  onset  of  a  PKP  wave  could  be  determined  to  ±0.5  sec 
and  compared  with  an  accurate  time  prediction  of  a  zero  depth  event,  the 
hypocenter  would  be  determined  to  ±  3  km.  (In  practice,  effects  of  local 
geology  cannot  be  completely  discounted.)  The  feasibility  of  this  is 
examined  briefly.  The  derived  curve  of  Figure  11  was  used  to  predict 
travel  times  for  surface  focus  events  recorded  at  Byrd  Station.  Appendix 
A  shows  the  results  of  subtracting  these  predicted  times  from  the  actual 
travel  times  without  depth  correction.  The  time  residual  can  be  used 
with  tables  from  Bolt  (1968)  to  recalculate  hypocenter  depth.  The  depths 
determined  from  our  residuals  are  then  compared  with  those  depths  given 
by  the  USC&GS.  All  those  events,  presumably  explosions,  listed  at  0  km 
by  the  USC&GS  also  show  times  corresponding  to  depths  within  5  km  of  the 
surface  based  on  the  Byrd  Station  arrival  times.  The  only  other  event 
with  such  a  shallow  depth  determination  is  an  event  at  Lake  Baikal 
(A  =  151°),  which  the  USC&GS  similarly  lists  as  4  km  deep. 

The  small  scatter  in  BC  arrival  times  allows  their  use  with  tele- 

seismic  P  in  least  squares  solutions  for  depth  determination.  The 

steeper  penetration  of  the  crust  and  upper  mantle  may  make  this  BC  curve 

more  valid  as  a  world  average  than  P  travel  time.  The  errors  shown  by 

Evernden  (1969)  in  depth  determination  from  average  P  travel  time  curves 

should  be  less  for  PKP. 

o 

Beyond  152  in  Appendix  A,  there  is  no  immediate  attempt  to  determine 
hypocenter  depth  from  arrival  times  because  of  differences  between  observed 
travel  time  and  published  values  (Bolt,  1968)  of  several  seconds.  BYA 
data  were  too  sparse  to  justify  fitting  a  curve  as  done  for  the  BC  branch. 
Arrivals  on  the  DF  branch  are  compared  to  the  Bolt  (1968)  curve.  The 
greater  scatter  observed  is  partly  due  to  interference  by  precursors. 
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For  A  >  160°,  on  larger  events  where  DF  and  AB  are  both  clearly 

observable,  there  Is  another  way  to  determine  hypocenter  depth  that  is 

Independent  of  the  USC&GS  epicenter  solutions.  Origin  times  are  still 

needed.  Figure  12  shows  the  recording  of  an  explosion  at  Novaya  Zemlya, 
o 

173  from  Byrd.  The  onset  of  both  DF  and  AB  arrivals  are  clear  on  all 
four  seismometers. 

At  A  ^  160°,  the  slope  of  the  AB  branch  is  many  times  as  steep  as 

that  of  the  DF  branch.  For  example,  the  interval  between  these  branches 

o 

is  increasing  at  4.0  sec/deg  at  173  .  If  the  AB-DF  interval  can  be 
picked  accurately  to  0.4  sec,  the  A  is  determined  to  0.1  deg.  Once 
this  confidence  can  be  gained  in  the  true  A,  the  depth  can  be  determined 
by  the  absolute  travel  time  and  predicted  travel  time  curves. 

Internal  Constitution  of  the  Earth 

Checks  were  made  for  azimuthal  dependence  of  residuals,  both  from 
teleseismic  P  and  PKP  arrivals;  no  dependence  was  found.  In  both  cases, 
the  majority  of  USC&GS-listed  events  occurred  in  a  small  range  of  azimuthal 
angles.  With  this  limitation,  a  small  variance  with  azimuth,  <  0.4  sec, 
would  not  have  been  detected. 

The  boundaries  of  the  outer  and  inner  cores  are  now  being  determined 

most  accurately  by  observation  of  PcP  and  PKIKP  waves.  Presently  used 

figures  for  their  mean  radii  are  3477  km  (Taggart  and  Engdahl,  1968)  and 

1216  km  (Bolt,  1968).  GH  arrivals,  which  could  confirm  the  extent  and 

velocity  of  an  intermediate  core,  are  not  clearly  enough  recorded  from 

the  small  events  used  in  this  study  for  quantitative  interpretation.  In 

agreement  with  Husebye  and  Toksoz  (1969),  we  find  GH  arrivals  cannot  be 

o 

observed  with  confidence  beyond  145  . 

The  observed  DF-AB  interval  in  A  >  160°  may  be  caused  by  deep 
variation  in  constitution  of  the  earth.  Bolt  (personal  communication) 
lacks  confidence  in  his  published  travel  times  as  a  world  norm  in  this 
range  due  to  systematic  anomalies  observed  in  the  region  of  his  test 
data  earthquakes  (Indian  Ocean).  Similarly,  of  course,  Byrd  Station  may 
exhibit  complimentary  anomalous  features.  An  alternate  hypothesis  put 
forth  is  that  the  core  of  the  earth  is  more  spherical  than  the  basic 
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ellipsoidal  figure  of  the  earth.  This  could  retard  AB  arrivals,  yet 
advance  DF  ones. 

The  arrival  of  energy  25  sec  before  the  AB  arrivals  at  A  =  173° 
might  be  interpreted  as  a  result  of  an  island  arc  crustal  plate  dipping 
into  the  mantle  at  Novaya  Zemlya.  In  view  of  similar  observations  from 
different  areas,  however,  we  cannot  strongly  support  this  conclusion. 

Figure  4,  shown  earlier,  plots  travel  time  from  well-recorded 
earthquakes  against  predicted  times.  Appendix  A  lists  the  events  used. 
Between  125  and  140°  there  were  many  more  arrivals  of  good  quality,  and 
only  the  clearest  onsets  were  used.  Outside  this  range,  the  number  of 
earthquakes  were  fewer,  and  consequently  all  events  which  appeared 
reliable  appear  on  the  figure.  Each  segment  of  the  travel  time  curve 
will  be  examined  separately. 

Arrivals  along  the  DE  segment  agree  closely  with  predicted  times 
o  o  o 

from  110  to  just  under  140  .  Near  140  ,  there  is  some  interference  by 
earlier  arrivals  of  increasing  amplitude  (see  Figure  8).  In  Appendix  A, 
these  early  arrivals  are  shown  without  depth  calculation.  This  inter¬ 
ference  may  hide  the  DE  onset  or  otherwise  lead  to  misinterpretation  of 
later  arrivals.  BYA  results  for  this  range  are  in  full  accord  with  the 
tabl es. 

Precursors  to  the  DE  arrivals  were  seen  on  large  earthquakes  from 

110  to  145°.  We  follow  Bolt's  nomenclature  and  call  these  GH.  Although 

these  precursors  tend  to  be  emergent  in  onset,  the  energy  arriving  before 

the  DE  arrival  increases  from  about  135°  up.  At  A  <  135,  the  GH  onsets 

are  seldom  picked  with  any  confidence.  At  A  >  142°,  there  is  some  question 

o 

in  distinguishing  GH  arrivals  from  other  energy.  Beyond  147  ,  where  time 

separation  would  again  permit  identification  of  GH,  the  arrivals  seen 

correspond  better  with  times  for  BC  (again  see  Figure  8).  BYA  data  can 

only  confirm  the  existence  of  precursors  to  the  DF  branch  and  give  only 

mildly  supporting  evidence  for  the  time  of  onset  of  this  branch. 

At  A  greater  than  145°,  the  EF  arrivals  were  recorded  from  very  few 

earthquakes.  Primarily,  there  simply  were  not  many  earthquakes  in  this 

region,  but  also  this  branch  no  longer  contains  the  majority  of  the 
o  .  o 

energy.  From  142  <  A  <  152  ,  the  BC  branch  dominates,  and  EF  arrivals 
show  up  as  precursors  only  on  some  of  the  larger  (mfa  5  5.5)  E.  Kazakh 
events.  At  A  >  152°,  the  separation  of  EF  from  other  arrivals  is  sufficient 
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to  avoid  interference;  however,  there  were  not  enough  large  earthquakes 
for  good  comparison  with  published  EF  values.  The  clearest  records  of 
EF  were  from  two  large  events  at  Novaya  Zemlya  (A  =  173°).  EF  arrivals 
from  these  were  2.1  and  2.6  sec  earlier  than  predicted  (Bolt,  1968). 
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SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

A  four-element ,  magnetic-tape-recording,  seismograph  array  was 
established  at  Byrd  Station,  Antarctica.  It  was  run  by  SRI  for  the 
calendar  year  1968  and  has  been  operated  since  by  USC&GS  personnel. 

Byrd  Station  was  found  to  be  a  fairly  quiet  site,  seismically, 
especially  within  the  band  of  PKP  frequencies.  Helicorder  records 
were  made  with  system  magnification  of  near  100  K.  The  limiting  factors 
were  the  station's  "cultural  noise''  and  low  frequency  background  noise, 
both  primarily  outside  the  signal  band.  The  detection  threshold  for  PKP 
events  from  unenhanced  BYA  records  is  shown  in  Figure  9.  No  significant 
correlation  was  found  between  noise  levels  and  local  storms.  Lack  of 
meteorological  data  prevented  correlation  with  oceanic  storms  to  the  north. 

Digital  tapes  were  made  of  142  events  for  examination  of  data 
enhancement  techniques.  Spectral  analysis  shows  the  PKP  signal  energy 
lies  primarily  between  0.4  and  2  Hz.  Digital  filtering  passing  this 
band  greatly  enhances  the  PKP  signals,  thereby  lowering  the  detection 
threshold.  Although  the  BYA  array  is  too  small  for  effective  velocity 
filtering,  the  straight  sum  of  signals  in  the  array  is  of  some  help  in 
observing  coherent  energy  arrivals.  The  data  now  on  hand  could  be 
enhanced  even  more  effectively  by  computerized  optimum  filtering. 
Digitizing  the  1969  data  would  roughly  double  the  data  sample  to  be 
treated. 

Work  with  records  from  BYA  for  the  first  two  years  showed  that  the 

BC  branch  of  PKP  is  sharp  enough  to  be  promising  for  depth  determination 

on  events  with  m.  >  4.8.  A  new  determination  of  detection  thresholds  as 
□ 

a  function  of  A  can  be  made  after  using  digital  techniques  to  enhance 
arrivals  on  more  records.  It  would  be  useful  to  work  out  an  analytical 
solution  for  the  possible  improvement  in  hypocenter  depth  determination 
by  adding  PKP  arrivals  to  the  least  squares  fit  of  teleseismic  P. 

The  EF  and  AB  branches  offer  possibilities  in  depth  determination 
o 

beyond  A  =  160  .  Arrays  in  South  Africa,  Malagasy,  and  Kerguelen  could 
be  used  to  monitor  western  U.S.  earthquakes  and  nuclear  shots  over  a 
range  of  A  to  establish  any  advantages  over  P  depth  determination.  Arrays 
should  be  large  enough  to  determine  dt/dA  as  an  aid  in  separating  the 
various  branches  near  BC. 
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Operating  arrays  covering  Asia  could  be  established  In  South  America 
and  New  Zealand,  even  if  only  the  BC  branch  should  have  a  clear  advantage 
over  P  in  depth  determination. 

Travel  times  computed  from  PKP  arrivals  observed  at  Byrd  show  general 
agreement  with  recently  published  tables  (Bolt,  1968).  Deviation  from 
the  most  recently  published  tables  include  the  following:  (1)  presence 
of  a  strong  BC  branch,  (2)  lack  of  confirmation  for  GH  branch  at 
A  >  145°,  (3)  greater  DF-AB  Interval  at  A  >  160°,  and  (4)  presence  of 
arrivals  between  the  DF  and  AB  arrivals  at  A  =  173°. 
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NOTES  FOR  APPENDIX  A 


0.  T.  origin  time  as  calculated  by  USC&GS 

*  Epicenter  determination  considered  by  USC&GS  as 

less  accurate  than  others. 

Lat.,  Long,  as  calculated  by  USC&GS. 

m,  body  wave  magnitude  as  calculated  bu  USC&GS. 

b 

No.  Sta.  number  of  stations  used  in  C&GS  computation. 

A°  geocentric  angle  between  seismometers  and  epicenter- 

average  used  when  pick  taken  from  more  than  one 
seismometer. 

Compression/dilation  -  first  motion  where  noted. 

T.T.  absolute  travel  time.  When  two  travel  times  are  listed, 

the  first  is  DF  and  the  second  is  either  BC  or  AB. 

ellipticity  correction-  from  Jeffreys  and  Bullen  1967. 

visual  appraisement. 

peak  amplitude  based  on  instrument  calibration  (two 
place  accuracy  only). 

X  indicates  digital  records  were  made  (by  Vela 
Digitizing  Center,  Alexandria,  Virginia). 

Residual  for  0  depth 

Predicted  time  minus  travel  time  where  predicted  time 
is  based  on  first  order  least  squares  fit  for  the  BC 
arrivals  (depth  corrected  using  A.  Qamar's  tables 
(Bolt,  1968)  for  Point  B  and  travel  times  have  been 
corrected  for  the  earths  ellipticity. 

Depth  from  Res. 

Using  A.  Qamar’s  tables  (Bolt,  1968)  for  Point  B 
C&GS  depth  -  33  km  is  used  where  C&GS  considers 
"n"  or  normal  depth. 

No.  of  seismometers 

the  number  of  seismometers  from  which  the  BYA  picks 
were  made  and  averages  taken.  The  1970  records  are 
generally  only  from  1;  the  helicorder  record.  Others, 

2,  3,  and  4  are  functions  of  playback  equipment  used. 
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